Non-Linear Vibrations in Nuclei 
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We have performed Time Dependent Hartree Fock (TDHF) calculations on the non linear response 
of nuclei. We have shown that quadrupole (and dipole) motion produces monopole (and quadrupole) 
oscillations in all atomic nuclei. We have shown that these findings can be interpreted as a large 
coupling between one and two phonon states leading to strong anharmonicities. 
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I. INTRODUCTION 

Fifty years ago, it was discovered that atomic nuclei 
may enter in resonance with electromagnetic fields [1]. 
This Giant Dipole Resonance (GDR) has been inter- 
preted as the vibration of neutrons against protons. Since 
then, other giant resonances (GR) have been predicted 
and observed, e.g. the Monopole GR (GMR), an alter- 
nation of compression and decompression of the nucleus, 
and the Quadrupole GR (GQR), an oscillation between 
a prolate and an oblate shape. The proof of the vibra- 
tional nature of the GR came only few years ago with the 
observation of the second vibrational quantum called the 
two-phonon state [2,3]. While many properties of these 
states plead in favor of an harmonic picture, striking ex- 
perimental observations such as an abnormally large ex- 
citation probability point to a strong coupling between 
the different phonon states [4,5]. This triggered a lot 
of theoretical investigations but only very weak anhar- 
monicities were found [6~11] creating an important crisis 
in our understanding of nuclear vibrations. 

Recently, it has been proposed that a strong anhar- 
monicity may come from large residual interaction lead- 
ing to the excitation of a GMR and a GQR on top of any 
state [12]. This was a surprise especially in the monopole 
case since it was generally believed that these couplings 
were small because of cancellation effects between var- 
ious diagrams [13]. It is shown in reference [12] that 
these cancellations between 3-particle 1-hole and 3-hole 
1-particle matrix elements were very limited. However, 
the approach of [12] even if it is fully microscopic do have 
some drawbacks. It is based on boson mapping meth- 
ods which may lead to violation of the Pauli principle 
and mixing with the spurious states [14]. Therefore, an 
independent confirmation of these important couplings 
leading to the excitation of a GMR and GQR on top of 
phonon states is crucial. 

In parallel looking to a completely different process, 
the excitation of a GDR in fusion reactions, we have 
shown [15] that, in a time dependent Hartree Fock 
(TDHF) approach, the dipole mode is non linearly cou- 
pled with other collective modes such as in particular the 
vibration of the density around a prolate shape. This 
work is also pointing in the direction of anharmonic vi- 
brations in nuclei but the particularities of the fusion 



dynamics and of the composite system does not allow 
to draw conclusions about the properties of the phonon 
built on the ground state. 

In this work, we present the first realistic TDHF cal- 
culation [16] of non-linear response to the collective vi- 
brations showing that, indeed, the one- to two- phonon 
coupling is a source of anharmonicities. We used the 
TDHF approach [17-19,16,20] which corresponds to an 
independent propagation of individual particles in the 
self-consistent mean field generated collectively. It does 
not incorporate the dissipation due to two-body interac- 
tion [21-23], but takes into account one body mechanisms 
such as Landau spreading and evaporation damping [24] . 
The quantal nature of the single particle dynamics is ex- 
plicitly preserved, which is crucial at low energy both 
because of shell effects and of the wave dynamics. In its 
small amplitude limit TDHF is equivalent to the Random 
Phase Approximation (RPA) which is the basic tool to 
understand the collective response of nuclei in terms of 
independent phonons. However, since the mean-field de- 
pends upon the actual excitation, TDHF is a non linear 
theory and hence contains couplings between collective 
modes. This point will be explicitly developed in the fol- 
lowing. In fact TDHF is optimized for the prediction of 
the average value of one body observables. Through non- 
linearities, it takes into account the effects of the residual 
interaction as soon as the considered phenomenon can 
be observed in the time evolution of a one body observ- 
able. Of course, the absence of terms explicitly taking 
into account the correlations is a limitation. In partic- 
ular, dampings and spreadings are neglected. As far as 
the time dependent approaches are concerned, it would 
be important to extend the present study to theories go- 
ing beyond the one body limit such as extended TDHF 
[23] which incorporate the effect of a " collision term" and 
also, through the fluctuations associated with the con- 
sidered dissipation, the coherent coupling with phonon 
plus particle hole excitations. Even more complete theo- 
ries such as the time dependent density matrix approach 
[10,25], which is known to reduce to the second RPA in 
its linearized version, would be an interesting extension 
of the present work. Finally, one should also try to apply 
the stochastic mean field approaches in particular in its 
version which have been proved to be potentially an exact 
solution of the many-body problem [26]. However, the 
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analysis presented in section 3 clearly show that the non 
linear response in TDHF contains the couplings between 
one and two phonon states coming from the 3-particle 
1-hole and 1-particlc 3-hole residual interaction. 

In section 2 we demonstrate first that coupling between 
one and two phonon states can be obtained through the 
evolution of average values of one-body obscrvablcs. In 
section 3, we show which part of the residual interaction 
is taken into account in a TDHF approach. In section 4 
we present results demonstrating the importance of the 
non-linear excitation of monopole and quadrupole modes 
on top of other collective vibrations. Finally we will con- 
clude in section 5. 

II. EFFECT OF COUPLINGS ON ONE-BODY 
OBSERVABLES 

To understand how this coupling can be extracted from 
the one-body dynamics let us consider the nonlinear cou- 
pling of a mode, | u > with the GR/i built on top of it 
leading to the two phonon state \ u/j, >. The Hamiltonian 
can be written 

H = Ho + V 

where Hq corresponds to the harmonic (RPA) part for 
which I z/ > and | vfi > are eigenstates with energies oJi, 
and ojiyf^ = uj^+ iVfj, while V is the residual interaction 
between phonons. For simplicity, let us introduce only 
the non-linear coupling =< v \ V \ vfi > which has 
been proven to be the most important one [12]. At the 
first order in £ = w^/o;^, this leads to the eigen states : 

\V >= {\ 1^ > -e\ ujj, >)/M 

and 

\vjl>= {£\v > -ir\v^ >)/Af 

where N"^ = 1 -|- A collective boost 

I ilj{t = 0) >= e-''^-'^" I - > 

inducing transitions between the ground state | — > and 
the collective state | u > with the amplitude Qi, =< — \ 

\ v >, leads to 

>^| - > -ik^q.e-"^-'\\ V > -e^'^^'e | TTJl >). (1) 

Then, < > (t) is simply given by 

<Q^> (t) ~ -2k^ql sm{io^t). (2) 

This shows that the linear response to collective boost 
induces oscillation of the collective moment at the col- 
lective frequency with an amplitude proportional to the 
transition probability q^. If we now compute the response 



to the operator which is associated with the excita- 
tion of the giant resonance fi, it is not zero because of the 
transitions between | V > and \i7Jl>. Using Eq. (1) and 
assuming < i^/i | | >= too, we get at the lowest 
order in k^, and v/tOf^ 

<Q^>{t)=2klqlq^^{cos{w^t)-l) (3) 

where the —1 term comes from higher order terms not 
explicitly written in Eq. (1). This demonstrates that the 
induced moment < > (t) is quadratic in the collective 
boost amplitude as expected from its non linear nature. 
Moreover, it oscillates at the frequency of the coupled 
mode with an amplitude proportional to the mixing 
coefficient v^/iUf^ and to the matrix clement q^. Finally, 
it should be noticed that < > (t) and < > (t) 
start in phase quadrature. 

III. LINEAR AND NON-LINEAR RESPONSE IN 
TDHF 

The TDHF approach is built to describe the average 
values of one-body observables. It propagates the evolu- 
tion of the one-body density matrix 

Pij =< ttj'ai > 

where is the operator creating a particle in the orbital 
\i >: 

idtp = [h, p] 

where h is the self-consistent mean-field Hamiltonian 
linked to the mean field energy E by hij = dE/dpji . 
We have used the code of ref [27] with SGII [28] and 
SLyAd [29] Skyrmc interactions. 

The RPA can be obtained by the linearization of the 
TDHF equation. Let us now go beyond the RPA by 
expanding the one-body density p up to the quadratic 
terms in the collective boost strength ki,: 
p = p(o)+k,pW+klp('), 

p^''^ the static HF groundstate, defines the occupied 
states (h) = J2h=i >< V'hl, the unoccupied 

states being the particle states (p). The condition on 
the one-body density p due to the independent particle 
approximation made when deriving the TDHF equation 
leads to the constrain = p. This imposes that p^^^ 
contains only p-h components which directly provide the 
particle-particle and hole-hole elements of the quadratic 
term of the one-body density 

J2) _ V- Jl) Jl) 
f^PP' ~ ^^ph fhp' 
h 

and 

(2) _ _ (1) (1) 

Phh' ^ / .PhyPyh'- 
P 
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A. TDHF and RPA 

The linear part of the TDHF equation leads to 

where 

is nothing but the RPA matrix acting only in the ph 
space. The RPA response after a boost e^'^'^"^" at time 
t = with Qi, = qvO^ + h.c. where 

is the RPA mode associated with the frequency cu^, is 

P^'^ = -g.p^ie-"^-* + h.c. . 

Using TrOyp+, = d^^', we get 

< Qv >= -2qlkvSmuj^t 

which corresponds to Eq. (2) and explain why the RPA 
provides a good approximation of and g^. 



where wc have explicitly factorizcd the transition proba- 
bility q^. Using TrO^p^, = (5^^/, the evolution of can 
be isolated: 

with 

< = TVO^VVq^,v^ = TrO^VVq^, 

and 

5v^ = TrO^^V/q2. 
The time independent parts of +Svtj. , v^, lead to 

so that 

2 

Comparing this result with the Eq. (3) shows that 
can be interpreted as the residual interaction exciting, in 
TDHF, the mode fi on top of the phonon v. 



B. Quadratic response and phonon coupling 

If we now compute the quadratic response in the 
TDHF approximation, we get for the p — h com- 
ponent of the one-body density p^^^ 

= A^p(2) +v^ + V^ + SV (4) 

where the three sources of non linearities are the p — h 
components of 

\V .W.W .(0)1 

ij k 

where £p = 1 and Sh = —1 and . It should be stressed 
that Eq. (4) is of course not the more general one but has 
the merit to be derived within the well defined framework 
of the time dependent mean field approximation. The 
p — h component of p^^^ can be expanded on the RPA 
basis 



C. Link with the residual interaction 

To illustrate this coupling, let us compute the contribu- 
tion coming from the forward amplitudes = O'^^^^ = 
Xpf^ then contains two terms involving the 3p — Ih 
and 3h — Ip residual interaction, 

'"ti = ^^ ^P'h";pp"X'^h^p"h"^p'h 

~ X] ^hh";h'p"XphXp»h"^ph' 

where we have introduced Viu^ji = dhji/dpki- Consider- 
ing the second term we get also two components 

= X] ^P'f>'";pp"^p"h"-^ph^p'h 

~ X] ^hh";h'p"^p"h"^ph'-^ph 

which are nothing but the exchange of jjL and v in v"^. 
These four terms correspond exactly to the X part of 
the phonon interaction V \ v) computed using boson 
mapping [12] except for a numerical factor which actually 
depends upon the mapping used. 

The X part of 5v, is E '""^Z'"' VW^P"h"- It 
results from the density dependence of the p — h inter- 
action defining the energy of the state v. In the simple 
case of a linear density dependence, it can be interpreted 
as a contribution of a three-body interaction inducing 
transitions from Ip — Ih to 2p — 2h. 
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This analysis clearly shows that the time dependence 
in TDHF takes into account the residual interaction. At 
the linear level, TDHF leads to the RPA. Going to the 
quadratic response, TDHF takes into account the one- to 
two-phonon coupling. The key point of the applicability 
of TDHF is the fact that the studied phenomenon can be 
deduced from the time dependence of the average value 
of a one- body observable. 



IV. RESULTS 

Let us now look at the TDHF results for the "^^Ca nu- 
cleus. We followed the monopole, quadrupole and dipole 
response for three initial conditions: 



A monopole boost using 
1 



Qo 



^{r^-<r^>{t = 0)). 



because of the spherical symmetry, a monopole 
boost can only trigger monopole modes. Therefore, 
we only observe < Qq > (t). 

A quadrupole boost generated by 

i 

The parity conservation forbids any dipole excita- 
tion when a quadrupole velocity field is applied to 
a spherical nucleus. Conversely, breathing modes 
(GMR) can be triggered by the quadrupole oscil- 
lation so that we do follow both the quadrupole 
< Q2 > {t) and the monopole < Qq > {t) re- 
sponses. 

An isovector dipole boost induced by 

Qd = Z/A ^Zn- N/A Zp. 



This excitation can be both coupled to the 
quadrupole and monopole oscillations so that we 
monitor the three moments, < Qo >(*)>< Q2 > 
{t) and <Qd> {t). 
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FIG. 1. Evolutions of the monopole, quadrupole and dipole 
moments (solid lines) and of the central density po (dashed 
lines) as a function of time for monopole (a), quadrupole (b) 
and dipole (c) excitations in "'"Ca. 

In figure 1, we observe that the collective boost induces 
oscillations of the associated moment as expected from 
the RPA (see Eq. (2)). They are only slightly damped in 
the GQR and GDR cases (fig. 1-b and 1-c respectively) 
while in the GMR case (fig. 1-a) beatings, characteristic 
of a Landau damping, are observed. This means that 
the dipole and quadrupole strengths are mostly concen- 
trated in a single resonance while the monopole one is 
fragmented. 

Plotting in figure 2 the amplitude of the first oscillation 
max a function of fc^ confirms the linearity of 
this response. Assuming that only one mode is excited 
which is a good approximation for the GDR and GQR 
(Eq. (2)) shows that the transition probability, ql, is < 
Qv >max /^k^- To get a deeper insight into the response 
we study the Fourier transform F{uj) of < Qu > {t)/k[y 
which is nothing but the RPA strength when the velocity 
field fcj/ is small enough to be in the linear regime. We 
see in figure 3 that the dipole and quadrupole modes are 
concentrated in a unique mode while the monopole is 
fragmented. However, the various peaks are in the same 
energy region so that they can be approximated by a 
single mode with a large Landau width. A detailed test 
of the equivalence between the linear regime of TDHF 
and the RPA response can be foimd in [24]. 

If we now turn to the non linearities, we can observe 
the moments which are different from the operator 
Qjy used for the exciting boost. We sec in figure 1 that, 
as expected from Eq. (3), this non linear response follows 
a {cos{u}ij,t) — 1) pattern oscillating with the frequency of 
the mode /i and not the one of the initially excited col- 
lective state i^. Moreover the amplitude of the first oscil- 
lation (fig. 2) is as expected quadratic in the excitation 
velocity /c^. In Fig. 1, one can see that large ampli- 
tude dipole (fig. 1-c) and quadrupole (fig. 1-b) motion 
induces variations of the central density po- Since the 
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central density can be modified only by monopole states 
this imposes that the large amplitude motion gets cou- 
pled with such breathing modes. In the same way a large 
amplitude dipolc oscillation induces a quadrupolc defor- 
mation of the nuclear potential and so gets coupled with 
the GQR. These observations lead to the conclusion that 
wc arc in the presence of a non-linear excitation of a gi- 
ant resonance ii on top of the collective motion v initially 
excited through the collective boost Qi,. 
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FIG. 2. Evolutions of the maximal oscillation amplitudes 
of < Qo > (solid lino), < Q2 > (dashed line) and < Qd > 
(dotted line) as a function of the intensities of the monopole 
(a), quadrupole (b), and dipole excitations (c). The horizon- 
tal lines represent the average number of excited phonons for 
each GR. 

The Fourier transform of < > (t) associated with 
the excitation of Qi, are also presented in figure 3. Let us 
first start with the quadrupole strength non-linearly ex- 
cited by a dipole boost. This is a clear indication that the 
observed state is indeed a GQR built on top of the GDR. 
This is what is expected from Eq. (3) where only the 
frequency LUf^ of the observed state appears. It should be 
notice that this frequency is different from one of the un- 
derlying dipole motion. The monopole case is more com- 
plex because of the presence of a strong Landau spreading 
and it seems that the strengths of the various monopole 
states depend upon the considered boost. This indicates 
that the coupling leading to the excitation of an addi- 
tional monopole state depends upon the collective mode 
initially excited. 
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FIG. 3. Monopole, quadrupole and dipole spectra (the 
Fourier transforms are in arbitrary units) for a monopole 
(left), quadrupole (middle) and dipole excitation (right). 

To estimate the magnitude of these non-linear cou- 
plings, we can first convert the amplitude of the induced 
oscillations into a phonon number using a coherent state 
picture. The horizontal lines in figure 2 represent the am- 
plitude of the oscillations associated with difiierent num- 
ber < rij/ >= 1, 2, 3, ... of excited phonons: 



<Q„> 



2 

max~ 



One can see that for a ki, which corresponds to the ex- 
citation of one phonon v the number of phonons fj, non- 
linearly excited is large. 

Assuming for each multipolarity a unique state |/U > 
non linearly excited one can use Eq. (3) to extract the 
residual interaction matrix element V/j, between > and 
l^lj. > from the amplitudes of the induced oscillations 

max 



Vn = <Qij. >max uJij./2klqlq^. 



(5) 



If the non-linear collective response is not concentrated 
in a unique state |/i > but corresponds to a set of 
states >} with uj^^ fa u>^, one can easily show that 
the extracted coefficient V/j, is related to the individual 
Uju^ =< i'\V\i'ij,i > by the weighted sum 

Vn ~ (6) 

where = qlJT,jql^ and g^, =< >. is 

in general higher than the individual w^^. For example, 
if the collective response is equally distributed into N 
states with identical coupling matrix elements ■ = t^^n 
then V/j, = VNv/j,. 
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TABLE I. Energies, transition probabilities q^, energy 

weighted sum-rules 5*1 and coupling coefficients of the GMR, 
GDR and GQR in the ^"Ca, ^°Zr and '^°'^Pb. and Si are 
expressed in fm? and MeV.frn^ for the GMR and GQR and 
in fm and MeVfrn^ for the GDR respectively. 



The results for the Ca, ^°Zr and '^'^^Pb arc presented 
in table 1. The lo^, are computed from the time to reach 
the first maximum of < Q,/ > (t). If the spreading of 
the observed mode is small, lu^ « m2/m,i which should 
be a little higher than mi /mo usually discussed. For the 
breathing mode, our results agree with the values ob- 
tained in [28] (mi /mo = 22.7, 19.5, 15.3MeV in "OCa, 
^°Zr and 208P6 respectively). These RPA results as well 
as our Wv are close to the averages which can be computed 
from the most collective states reported in ref. [12]. For 
the '*°Ca, the m2/mi for the monopole, quadrupole and 
dipole states of ref. [12] are respectively 21.2, 16.9 and 
18.5 MeV for 84%, 85% and 66% of the corresponding 
energy weighted sum rules (EWSR). For 208Pb these val- 
ues are 14.1, 11.1 and 13.6 MeV for 89%, 91% and 80% 
of the EWSR respectively. The relative sign of and 
is given by the early evolution of the moments described 
by Eq. (5) and shown in figure 1. They appear to be 
all negative in agreement with ref. [12]. The couplings 
are large of the order of few MeV. From the quan- 
titative point of view, the non-linear coupling extracted 
from TDHF appears to be 50% larger than the one re- 
ported in reference [12]. This is a reasonable agreement 
since TDHF result is a weighted sum of the individual 
couplings as shown in Eq. 6. Summing the contribu- 
tions of the different collective states considered in ref. 
[12] reduces the difference between the reported values. 
However, the phonon basis studied in ref. [12] being in- 
complete it is expected that the TDHF results remains 
higher. It should be also noticed that some difference can 
remain due to the approximations involved in the differ- 
ent approaches as discussed in the quadratic response 
analysis. 

In table 1 one can also see that the larger the nucleus 
the smaller the coupling. This is in agreement with the 
fact that these couplings are mediated by the surface. 
To control the robustness of our conclusion we have per- 
formed a series of calculations using a different Skyrme 
force, the recent SLy4d parametrization. For a ^°Ca, this 
leads to a coupling exciting the GMR on top of the GDR 
of -4.01 MeV and of -4.36 MeV on top of the GQR. 
The quadrupole response during a dipole oscillation leads 
to a residual interaction of —3.98 MeV. Those results are 
very close to the one reported in table 1. 



V. CONCLUSIONS 



In conclusion, we have shown with TDHF calculations 
that a non-linear excitation of monopole and quadrupole 
should occur on top of any collective motion in nuclei. 
These couplings can be interpreted in terms of a large 
residual interaction which couples one-phonon and two- 
phonon states. These results show that large anhar- 
monicities should be expected in the collective motions 
in nuclei. 
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